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O ; ABSTRACT 

\^ • Combining cosmological SPH simulations with a previously developed Lya produc- 

^ (— I tion/transmission model and the Early Reionization Model (ERM, reionization ends 

O ^ at rcdshift z ~ 7), we obtain Lya and UV Luminosity Functions (LFs) for Lyman 

Alpha Emitters (LAEs) at 5.7 ^ z ^ 7.6. Matching model results to observations at 
z ~ 5.7 requires escape fractions of Lya, f a = 0.3, and UV (non-ionizing) contin- 
uum photons, f c = 0.22, corresponding to a color excess, E(B — V) = 0.15. We find 
that (i) f c increases towards higher redshifts, due the decreasing mean dust content 
of galaxies, (ii) the evolution of f a /fc hints at the dust content of the ISM becoming 
progressively inhomogeneous/clumped with decreasing redshift. Using the model as- 
s ! ' sumptions, clustering of sources has little effect on the Lya LF for a cosmic hydrogen 

neutral fraction \hi ^ 10~ 4 , a value attained at z ^ 6.6 in the ERM. However, during 
the initial reionization phases (z XL 7) the clustering photoionization boost becomes 
important. We quantify the physical properties of observed LAEs and their redshift 
n . evolution, for which we give handy analytical fitting functions. Halo (stellar) masses 

are in the range 10.0 < \ogM h < 11.8 (8.1 < log A/* < 10.4) with M h oc Mf 64 . The 
star formation rates are M* = 3 — 120 M Q yr" 1 , mass- weighted mean ages are t» > 20 
Myr at all redshifts, while the mean stellar metallicity increases from Z = 0.12 Z Q at 
\ z ~ 7.6 to Z = 0.22 Z Q at z ~ 5.7; both t* and Z positively correlate with stellar 

mass. The brightest LAEs are all characterized by large M* and intermediate ages 
(« 200 Myr), while objects in the faint end of the Lya LF show large age and star 
formation rate spreads. With no more free parameters, the Spectral Energy Distribu- 
jJJ ■ tions of three LAE at z ~ 5.7 observed by Lai et al. (2007) are well reproduced by an 

intermediate age (182 — 220 Myr) stellar population and the above E(B — V) value. 
The model uncertainties, mostly related to the simplified treatment of dust and to the 
possible effects related to gas outflow/infall, are discussed along with their impact on 
the results. 

Key words: mcthods:numerical - galaxies:high redshift - luminosity function - in- 
tergalactic medium - cosmology: theory 
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1 INTRODUCTION 

Over the past few years, Lyman Alpha Emitters (LAEs) 
have rapidly been gaining popularity as probes of cosmic 
reionization for two primary reasons. Firstly, specific signa- 
tures like the strength, width of the Lya line (1216 A) and 
the continuum break bluewards of it make the detection of 
LAEs unambiguous. Secondly, since Lya photons are highly 
sensitive to the presence of neutral hydrogen, their attenu- 
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ation can be used to put constraints on the ionization state 
of the IGM. 

Since these are also amongst the earliest galaxies to 
have formed, they represent superb probes of the proper- 
ties and evolution of early galaxy populations. This is im- 
portant since mechanical, chemical and radiative feedback 
from these galaxies will determine the properties of galaxies 
formed later on. Obtaining the star formation rates (SFR), 
metallicity, initial mass function (IMF) of these galaxies are 
hence imperative in understanding galaxy evolution. 

Advances in instrument sensitivity have enabled ob- 
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servers to push the observable frontier to increasingly high 
redshifts. There are now hundreds of confirmed LAEs at 
z ~ 2.25 (Nilsson et al. 2008), z ~ 3 (Cowie & Hu 1998; 
Steidel et al. 2000; Matsuda et al. 2005; Venemans et al. 
2007; Ouchi et al. 2008), z ~ 4.5 (Finkelstein et al. 2007), 
z ~ 5.7 (Malhotra et al. 2005; Shimasaku et al. 2006) and 
z ~ 6.6 (Taniguchi et al. 2005; Kashikawa et al. 2006). 

While the data accumulated on LAEs shows no evolu- 
tion in the apparent Lya luminosity function (LF) between 
z = 3.1 - 5.7 (Ouchi et al. 2008), the LF changes appreciably 
between z = 5.7 and 6.6 (Kashikawa et al. 2006) with L* 
at z = 6.6 being about 50% of the value at z = 5.7. Sur- 
prisingly, however, the UV LF does not show any evolution 
between these same redshifts. Kashikawa et al. (2006), have 
proposed this evolution in the Lya LF to be indicative of a 
sudden change in the ionization state of the Universe. How- 
ever, the problem of why reionization would affect the high 
luminosity end of the LF rather than the the low, faint end, 
as expected, remains. 

A range of theoretical models, both semi-analytic (Di- 
jkstra et al. 2007a,b; Kobayashi et al. 2007,2009 ; Dayal, 
Ferrara & Gallerani 2008 [DFG08]) and those involving cos- 
mological simulations (McQuinn et al. 2007; Nagamine et al. 
2008) have been built to explain the observations and use 
LAEs as probes of reionization. Using a search in a param- 
eter space comprised by the SFR efficiency and Lya trans- 
mission, Dijkstra et al. (2007b), showed that the evolution 
in the LF between z = 5.7 and 6.6 could be explained solely 
by an evolution of the underlying mass function. Kobayashi 
et al. (2007, 2009), used a semi-analytic galaxy formation 
model, including the effects of dust and outflows on the es- 
cape fraction of Lya photons to explain the same. Nagamine 
et al. (2008), instead, have used the duty cycle scenario (the 
fraction of LAEs turned on at a certain time) to explain the 
UV LF at 2 = 3 and 6. McQuinn et al. (2007), have shown 
that reionization increases the measured clustering of LAEs 
and hence, these objects are very useful to probe the epoch 
of reionization. 

In DFG08 we showed that the Lya and UV LFs can 
be explained by an evolution of the underlying dark matter 
halo mass function and dust attenuation that increases with 
decreasing redshift. We found that the ERM (Early Reion- 
ization Model), in which reionization ends at z = 7, explains 
the data more consistently as compared to the LRM (Late 
Reionization Model, reionization ends at z = 6). We thus 
found that reionization does not play any role in shaping the 
Lya LF at z < 6.6. This is consistent with the results ob- 
tained by McQuinn et al. (2007), who find that the Universe 
must be highly ionized at z ~ 6.6. In addition to the LFs, 
our model also reproduced the weighted skewness measure- 
ments (Kashikawa et al. 2006), the equivalent width (EW) 
at z — 4.5 (Dawson et al. 2007) and predicted the SFR den- 
sity of LAEs. However, both the average and distribution of 
EWs from our model at z = 5.7 were much lower than those 
observed by Shimasaku et al. (2006). 

In spite of these different approaches, there are still a 
number of missing ingredients, the most important of these 
being the calculation of the intrinsic Lya luminosity, UV lu- 
minosity and SED (spectral energy distribution) as a func- 
tion of the SFR, age, metallicity and IMF of the galaxy un- 
der consideration. Other important effects include the boost 
in the ionization rate imparted by galaxy clustering, its ef- 



fects on the visibility of galaxies of different masses and the 
effects of inflows/outflows on the Lya luminosity. Up to now, 
only a general value of this boost has been used (Dijkstra et 
al. 2007a); however, the effects of clustering on the visibil- 
ity of LAEs as a function of galaxy properties have largely 
remained unexplored. 

In this paper, we use state-of-the-art cosmological SPH 
simulations to fix the SFR, age, metallicity for each galaxy 
to obtain the intrinsic Lya luminosity, UV luminosity and 
the SED. We again use the ERM (Gallerani et al. 2007), 
which accounts for all available data beyond LAEs, including 
Lya/Ly/3 Gunn- Peterson opacity, electron scattering optical 
depths, Lyman limit systems, cosmic SFR histories and the 
number density of high redshift sources. Using the above 
ingredients, we obtain the Lya and UV LFs and we are 
able to quantify the importance/effect of clustering on Lya 
luminosity transmission and its contribution to shaping the 
Lya LF. By doing so, we gain insight on the nature of LAEs 
and put constraints on their elusive physical properties. 



2 SIMULATIONS 

The simulation analyzed in this paper has been carried out 
using the TreePM-SPH code GADGET-2 (Springel 2005), 
with the implementation of chemodynamics as described by 
Tornatore et al. (2007). It is part of a larger set of cos- 
mological runs, which are presented and discussed in de- 
tail elsewhere (Tornatore et al. 2009, in preparation). The 
adopted cosmological model corresponds to the ACDM Uni- 
verse with fi m = 0.26, «a = 0.74, fi b = 0.0413, n s = 0.95, 
Ho — 73 km s _1 Mpc -1 and erg = 0.8, thus consistent 
with the 5-year analysis of the WMAP data (Komatsu et al. 
2009). The periodic simulation box has a comoving size of 
75/i _1 Mpc and contains 512 3 Dark Matter particles and ini- 
tially, an equal number of gas particles. As such, the masses 
of the DM and gas particles are hidm — 1.7 x lO 8 /i _1 M 
and m gas — 4.1 x 1O 7 /i _1 M , respectively. The Plummer- 
equivalent softening length for the gravitational force is set 
to epi = 2.5 /t _1 kpc, kept fixed in physical units from z — 2 
to z = 0, while being epi = 7.5 h~ kpc in comoving units 
at higher redshift. The value of the softening parameter of 
the SPH kernel for the computation of hydrodynamic forces 
is allowed to drop at most to half of the the gravitational 
softening. 

The run assumes a metallicity-dependent radiative cool- 
ing (Sutherland & Dopita 1993) and a uniform redshift- 
dependent Ultra Violet Background (UVB) produced by 
quasars and galaxies as given by Haardt & Madau (1996). 
The code also includes an effective model to describe star 
formation from a multi-phase interstellar medium (ISM) and 
a prescription for galactic winds triggered by supernova (SN) 
explosions (see Springel & Hernquist 2003 for a detailed de- 
scription). Galactic winds are assumed to have a fixed ve- 
locity of 500 km s _1 , with a mass upload rate equal to twice 
the local star formation rate. The code includes the descrip- 
tion of chemical enrichment given in Tornatore et al. (2007) . 
Metals are produced by SNII, SNIa and intermediate and 
low-mass stars in the asymptotic giant branch (AGB). We 
assume SNII arise from stars having mass above 8M©. As 
for SNIa, we assume their progenitors to be binary systems, 
whose total mass lies in the range (3-16) M Q . The relative 
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number of stars of different mass is computed for this simu- 
lation by assuming the Salpeter (1955) IMF between 1 and 
100 Mq. Metals and energy are released by stars of differ- 
ent masses by properly accounting for mass-dependent life- 
times. In this work we assume the lifetime function proposed 
by Padovani & Matteucci (1993). We adopt the metallicity- 
dependent stellar yields from Woosley & Weaver (1995) and 
the yields for AGB and SNIa from van den Hoek & Groe- 
newegen (1997). 

As for the identification of galaxies, they are recognized 
as gravitationally bound groups of star particles. For each 
analyzed snapshot we first run a standard friends-of-friends 
(FOF) algorithm with a linking length of 0.2 in units of the 
mean particle separation. Each FOF group is then decom- 
posed into a set of disjoint substructures, which are identi- 
fied by the SUBFIND algorithm (Springel et al. 2001) as lo- 
cally overdense regions in the density field of the background 
main halo. After performing a gravitational unbinding pro- 
cedure, only sub-halos with at least 20 bound particles are 
considered to be genuine structures (see Saro et al. 2006, for 
further details). For each "bona-fide" galaxy, we compute 
the mass- weighted ag^E the total halo/stellar/gas mass, the 
SFR, the mass weighted gas/stellar metallicity, the mass- 
weighted gas temperature and the half mass radius of the 
dark matter halo. 

We compute the Lya emission and spectral properties 
for all the structures identified as galaxies in the simulation 
boxes at the redshifts of interest (z ~ 5.7,6.6,7.6). Obvi- 
ously, not all these galaxies will be necessarily classified as 
LAEs. 



3 PHYSICS OF LYMAN ALPHA EMISSION 

We summarize here the main features of our model, the de- 
tails of which can be found in DFG08. Star formation in 
galaxies gives rise to continuum and H i ionizing photons, of 
which the the latter ionize the Inter Stellar Medium (ISM). 
Due to the high density of the ISM, recombinations take 
place on a short timescale and this gives rise to a Lya emis- 
sion line. 

The intrinsic UV continuum and Lya luminosity de- 
pend on the galaxy properties including the IMF, SFR 
(M„), stellar metallicity (Z) and age (£,); all these quan- 
tities are taken from the simulation outputs, as discussed 
above. STARBURST99, (Leitherer et al. 1999), a popula- 
tion synthesis code, is used to obtain the intrinsic values of 
the continuum luminosity as well as the rate of H i ionizing 
photons produced using the Salpeter IMF, appropriate SFR, 
age and metallicity for each simulated galaxy. From our 
model, for a galaxy with U = 200 Myr, Z = O.2Z , 
M* = 1M© yr _1 , the rate of production of H 1 ionizing pho- 
tons, Q — 10 53,47 s _1 and the corresponding intrinsic Lya 
luminosity is L™' = 3.25 x 10 42 ergs _1 . This is very consis- 
tent with the value of 3.3 x 10 42 ergs _1 shown in Tab. 4 of 
Schaerer (2003) for similar values of the age, metallicity and 



IMF. For the same galaxy, the intrinsic continuum luminos- 
ity is Ll nt = 3.5 x lO^ergs^A" 1 , which yields an EW of 
about 93 A. 

The Lya line profile is affected by the rotation velocity 
of the galaxy (v c ), unlike the continuum. For quiescent star 
formation, for realistic halo and disc properties, the galaxy 
rotation velocity can have values between 1-2 times the halo 
rotation velocity (Mo, Mao & White 1998; Cole et al. 2000). 
We use a value of 1.5 in this work. We calculate the velocity 
of the halo, Vh, assuming that the collapsed region has an 
overdensity of roughly 200 times the mean cosmic density 
contained in a radius r2oo- Then, Vh, the velocity at r^oo is 
expressed as 
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where M^ is the dark matter halo mass and Q m (z), H(z) 
are the matter density and Hubble parameters, respectively, 
at the redshift of the emitter. 

Two processes determine the Lya luminosity which 
emerges from the galaxy. First, only a fraction of the 
H 1 ionizing photons ionize the ISM, contributing to the Lya 
luminosity while the rest (/esc) escape the galaxy and ionize 
the Inter Galactic Medium (IGM) surrounding it. Second, 
only a fraction, f a , of the Lya photons produced escape the 
galaxy, unabsorbed by dust in the ISM. 

The Doppler broadened Lya luminosity profile that 
emerges from the galaxy is expressed as 



L e a m (v) = -Qhv a {l - f e3c )f a —=r- — 
3 ^/ivAud 



exp 



(2) 



where Q is the H 1 ionizing photon production rate for the 
intrinsic galaxy properties, h is Planck's constant, v a is the 
frequency of Lya photons, Aud ~ [v c /c]v a , c is the speed 
of light and the factor two-thirds arises assuming Case B 
recombination (Osterbrock 1989). 

The continuum band (1250-1550 A, centered at 1375 
A) is chosen such that it is unaffected by attenuation due 
to H 1 and it is only attenuated by dust in the ISM. The 
fraction, f c (f a ) of continuum (Lya) photons which escape 
the galaxy depends on the total amount of dust in the ISM 
as well as its distribution (homogeneous or clumped). The 
continuum luminosity emerging from the galaxy is expressed 

„„ r em _ r r int 

Lya photons suffer further attenuation as they travel 
through the IGM due to their large H 1 optical depth; even 
small amounts of neutral hydrogen in the IGM can attenu- 
ate the Lya luminosity by large amounts. To calculate this 
attenuation, we use the prescription detailed in Sec 2.4 of 
DFG08. We use the value of the UVB photoionization rate 
(r_g) from Gallerani et al. (2007) for the Early Reioniza- 
tion Model (ERM), wherein reionization ends at z ~ 7, to 
calculate the mean neutral hydrogen fraction in the igm3, 
■% HI = riHi/nH- In principle, our simulation could be used 
to account for inhomogeneities in the IGM along differ- 
ent lines of sight. However, this work is much beyond the 
scope of this paper and so, we assume a homogeneous and 



1 This method tends to slightly bias the age towards larger val- 
ues. On the other hand numerical resolution limits the ability to 
resolve the smallest halos harboring the oldest stars. 



2 This is slightly inconsistent with the UVB used in the simula- 
tions. However, this produces negligible differences in the proper- 
ties of the LAEs. 
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isotropic IGM density field. According to these calculations, 
at z ~ (5.7, 6.6, 7.6), X m = (6 X 10" 5 , 2.3 x 10 -4 , 0.16). For 
the given value of xhi, we calculate the volume (Vr) of the 
so-called Stromgren sphere built by the H i ionizing photons 
that escape the galaxy without ionizing the ISM and instead 
ionize the IGM (Shapiro & Giroux, 1987; Madau, Haardt & 
Rees, 1999) as: 



dVi 
dt 



3H(z)Vi 
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XHinn t rec 
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where jih is the hydrogen number density at the redshift 
considered and t rec is the volume averaged recombination 
time in the IGM. 

Estimating f esc proves to be somewhat of a challenge. 
While empirical estimates from normal local and high red- 
shift galaxies have yielded modest values of a few percent, 
theoretical studies have been largely inconclusive. Recently, 
using an N-body+ hydrodynamic simulation, Gnedin et al. 
(2008), have found an average value, / esc ~ 0.02, for halo 
masses between 1O 1C) - 12 M over redshifts between 3 and 5. 
We use this value as a reasonable best estimate of / esc for 
all the galaxies in the simulation boxes at z ~ 5.7, 6.6, 7.6 to 
calculate the size of their Stromgren Sphere. We could also 
calculate f eac by coupling our simulation with one that deals 
with radiative transfer. However, such a detailed analysis is 
deferred to future works. 

Inside the Stromgren sphere, the ionization rate has two 
contributions : (a) a constant value, Tb, determined by the 
UVB photoionization rate and (b) a radius dependent value, 
Tl , determined by the luminosity that emerges out from the 
galaxy. The total photoionization rate (Tbl) at any distance 
r from the galaxy can be expressed as 

r BL (r) = r s + r L (r), 
r BL (r) = t b + 



where Z^ m = L^'/esc [ergs -1 Hz -1 ], is the specific ioniz- 
ing luminosity emerging from the emitter, vl is the fre- 
quency corresponding to the Lyman limit wavelength (912 
A) and ox is the hydrogen photoionization cross-section. 
Since the IGM is in local photoionization equilibrium, we 
use ionization-recombination balance to compute the value 
of xhi at each point within the Stromgren Sphere. At the 
edge of this region, we force xhi to attain the constant 
global value, i.e., Tbl ~ Tb. 

If z em and z bs are the redshifts of the emitter and ob- 
server respectively, we calculate the total optical depth (r a ) 
to the Lya photons along the line of sight (LOS) as 
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where dl/dz = c/[H(z)(l + z)]. Here, a is the total ab- 
sorption cross-section and <f> is the Voigt profile. We use 
do = ire 2 f / (m e c) , where e, m e are the electron charge and 
mass respectively and / is the oscillator strength (0.4162). 

For regions of low H i density, the natural line broad- 
ening is not very important and the Voigt profile can be 
approximated by the Gaussian core: 




Figure 1. A schematic diagram to elucidate the effects of cluster- 
ing on galaxy visibility (see text). The circles indicate the bound- 
aries of individual H II regions. 



b * 



In eq. [6] Vi is used since a photon of initial frequency v 
has a frequency Vi — v[(l + Zi)/(1 + Zem)] at a redshift Zi 
along the LOS. The Doppler width is expressed as b/X a , 
where b — yj2kT /m,H is the Doppler width parameter, ran 
is the hydrogen mass, k is the Boltzmann constant and 
T = 10 4 if is the IGM temperature (Santos 2004; Schaye 
et al. 2000; Bolton & Haehnelt 2007). In more dense regions 
the Lorentzian damping wing of the Voigt profile becomes 
important. According to Peebles (1993), this can be approx- 
imated as 



(7) 



where A = 6.25 x 10 8 s -1 is the decay constant for the 
Lya resonance. Although computationally more expensive 
than the above approximations, using the Voigt profile to 
compute the absorption cross-section gives precise results, 
and therefore we have implemented it in our code to obtain 
all the results presented below. 

The observed Lya luminosity (La) can be expressed as 
L a — L^Ta where a fraction T a = e~ Ta of the Lya lu- 
minosity emerging from the galaxy is transmitted through 
the IGM. Since the continuum is unaffected by transmis- 
sion through the IGM, the observed continuum luminosity 
is calculated to be L c = L E r m . 



4 EFFECTS OF CLUSTERING 

One of the most important drawbacks of modelling LAEs 
semi-analytically is that only an average value of the UVB 
photoionization rate, r_g, is usually assumed, making any 
inference on the additional contribution, Tg, due to galaxy 
clustering not possible. As McQuinn et al. (2007), have 
shown, the observed clustering of LAEs is increased by trans- 
mission through a patchy IGM as expected during reioniza- 
tion. In addition to being an important tool to probe reion- 
ization, clustering will also leave imprints on the Lya LF. 

Thanks to our cosmological SPH simulations, we are 
able to include the effect of clustering on the Lya LF and 
transmission. Clustering of galaxies implies that the mean 
separation could become smaller than the typical size of 
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their ionized regions. Therefore, more than one galaxy af- 
fects the size of the H n region and the neutral hydrogen 
profile within it, thereby increasing the visibility of indi- 
vidual galaxies. Consider for example the case of a luminous 
galaxy, A, and a fainter one, B (see Fig. 1). Each one of them 
will carve an H n region in the IGM, whose size depends on 
their luminosity, as shown in the Figure. Depending on the 
IGM hydrogen neutral fraction, none, one or both galaxies 
can be seen individually. This is because, to a first approx- 
imation, the spatial scale imposed by the Gunn-Peterson 
damping wing on the size of the H n region corresponds 
to a redshift separation of Az ~ 0.01, i.e. about 200 kpc 
(physical) at z = 10 (Miralda-Escude 1998). However, when 
the mean separation between galaxies is smaller than the 
size of the smallest H n region (i.e. galaxies are consider- 
ably clustered) then it can occur that both galaxies can be 
observed as LAEs due to the increased H n region size and 
decreased H i fraction within it boosting the transmissivity. 
The importance of such clustering effect on LAE visibility 
is stronger in the initial reionization phases, when the gas is 
almost neutral, implying that only clustered galaxies would 
be visible under those circumstances. 

To quantify this effect from simulations, we use a post- 
processing technique in which we start by calculating the 
size of the ionized region around each galaxy as explained in 
eq. ©. Any two galaxies that are separated by a distance 
smaller than either of their Stromgren radii are then treated 
as a local enhancement of the photoionization rate for the 
other, i.e., for galaxy A, we add an extra contribution, To, 
evaluated as the local value of the photoionization rate of 
galaxy B at the position of galaxy A. The same procedure 
is followed when considering the transmissivity of galaxy B. 
More generally, the total ionization rate, Ft, seen at the po- 
sition of galaxy j, whose separation from N other galaxies is 
smaller than either of the Stromgren radii can be expressed 
as 

A poo j em / \ 3 

i* jM -r„ W + £ / (» 

where the second term on the right hand side represents To- 
Further, LJ™ is the ionizing luminosity emerging from the 
i th galaxy and Tij is the radial distance between galaxies i 
and j. As a caveat, we point out that we are calculating clus- 
tering effects assuming that the sizes of the H n regions cor- 
respond to sources embedded in the IGM at the mean ioniza- 
tion fraction given by the ERM. However, prior to complete 
overlap, the ionization field is very patchy and galaxies are 
more likely to be immersed in either an almost neutral or 
a highly ionized region. Because the average \BI = 0.16, it 
means that 1 — \hi = 84% of the volume is substantially 
ionized at z = 7.6. Hence, in general, our method provides 
a lower limit to the number of detectable LAEs. Several au- 
thors (Zahn et al. 2007; Mesinger & Furlanetto 2007; Geil 
& Wyithe 2008) have presented schemes that avoid detailed 
radiative transfer (RT) calculations and still provide ioniza- 
tion schemes in good agreement with simulations. A precise 
calculation of clustering effects can also be done by properly 
following radiative transfer in detail. We defer these inclu- 
sions to our model to further work. 

To summarize, the boost in the ionization background 
imparted by clustering is important for all galaxies when the 



IGM is close to neutral. The importance of this effect de- 
creases as reionization proceeds, however, it decreases faster 
for more luminous galaxies. 



5 RESULTS 

We are now ready to compare model results with obser- 
vations. In particular, we compare the calculated UV and 
Lya LFs to the data obtained by Shimasaku et al. (2006) 
(z ~ 5.7) and Kashikawa et al. (2006) (z ~ 6.6). We also pre- 
dict the LFs for z ~ 7.6 although data are not yet available 
at this redshift. In addition, we also explore the effects of 
clustering on Lya transmission and present synthetic spec- 
tral energy distribution (SEDs) to match with Lai et al. 
(2007) observations at z ~ 5.7. Finally, we will discuss the 
physical properties such as age, metallicity, halo/stellar mass 
and SFR for all objects identified as LAEs from this work. 

5.1 Lya LF and the effects of clustering 

As mentioned before, we use the procedure explained in 
Sec. [3] to calculate the observed Lya luminosity for each 
of the galaxies in the simulation boxes. Galaxies with (a) 
observed Lya luminosity in the currently observable range, 
L a ^ 10 42 ' 2 ergs _1 and (b) value of the observed EW, 
EW > 20 A are then identified as LAEs, which are used 
to build the cumulative Lya LF to compare to observations. 
The number of objects identified as LAEs from our simu- 
lations are (1696, 929, 136) at z ~ (5.7,6.6,7.6); since the 
definition of LAE is an operational one based on the ob- 
served Lya luminosity, the number depends on the adopted 
values of f a and T a . The values given above are for the best 
fit parameters shown in Tab. 1 and Fig. [2] 

We find that, independent of clustering, to match the 
data at both z ~ 5.7 and 6.6, only a certain fraction 
(fa = 0.3) of the Lya luminosity must emerge out of the 
galaxy, unabsorbed by dust within the ISM. As shown in the 
uppermost panel of Fig. within error bars, the theoretical 
LF nicely matches the data at z ~ 5.7. 

With the same value of f a — 0.3 at z ~ 6.6, the theoret- 
ical LF lies close to the upper bound of the data as shown in 
the central panel. Since the dust attenuation of Lya photons 
does not evolve between z ~ 6.6 and 5.7, it is reasonable to 
assume that the same value also holds at z ~ 7.6 (the cosmic 
time between z ~ 7.6 and 6.6 is only about 0.1 Gyr, which 
might be too small for any significant dust evolution) and 
then use it to predict the LF (bottom panel of Fig. ■ Fig. [2] 
provides interesting information on the effects of clustering 
as well: from the overlapping of the LFs including/excluding 
clustering at z ~ 5.7 and 6.6, we conclude that these are neg- 
ligible at these epochs. Clustering of sources, however, plays 
a key role at z ~ 7.6. If clustering effects are neglected, very 
few objects (about 30) would be luminous enough to be 
detectable in current Lya surveys; instead, the luminosity 
boost due to clustered sources leads to about 136 objects to 
become visible in our simulation volume. We reiterate that, 
because we do not follow the radiative transfer of ionizing 
radiation through the IGM, the results presented in Fig. 2 
(and in Fig. 4) at z = 7.6 must be seen as lower limits to 
the actual number density of LAEs. 

We pause to discuss an issue concerning the LF at 
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Table 1. Best fit parameter values to match the Lya and UV LFs in the ERM, including clustering effects. For each redshift (col. 1), 
we report the UVB photoionization rate (col. 2), the fraction of neutral hydrogen corresponding to the UVB photoionization rate (col. 
3), the average value of the boost parameter (col. 4), the average transmission of the Lya luminosity (col. 5), the value of the the escape 
fraction of Lya photons (col. 6), the escape fraction of UV continuum photons (col. 7) and the color excess calculated using the supernova 
dust extinction curve (col. 8). 
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Figure 2. Cumulative Lya LF for the ERM. The panels are for 
z ~ 5.7, 6.6 and 7.6 from top to bottom. Points represent the 
data at two different redshifts: z ~ 5.7 (Shimasaku et al. 2006) 
(circles), z ~ 6.6 (Kashikawa et al. 2006) with downward (up- 
ward) triangles showing the upper (lower) limits. Solid (dashed) 
lines in each panel refer to model predictions at z ~ 5.7, 6.6, 7.6 
for the parameter values in Tab. [l] including (excluding) cluster- 
ing effects. Shaded regions in all panels show poissonian errors. 
Though curves with/ without clustering effects are indistinguish- 
able at z = 5.6, 6.6; the difference is appreciable at z = 7.6. Be- 
cause radiative transfer is not accounted for in the simulations, at 
2 = 7.6 our results represent lower limits to the observed objects. 



2 ~ 6.6. The model predictions for the faint end of the LF 
lie between the upper and lower limits set by the observa- 
tional data. Although this could be regarded as a success 
of the model, it is safe to discuss if physics not included in 
our model or biases in the data might spoil this agreement. 
While the lower bound of the data is made up of spectro- 
scopically confirmed LAEs, the upper bound is a photomet- 
ric sample composed of all galaxies identified as LAE candi- 
dates at this redshift. Spectroscopic analysis of LAE candi- 
dates on the upper bound of the LF could then possibly rule 
out a number of them as being low-z interlopers and con- 
taminants, thus bringing the upper limit in agreement with 



our curve. This would require that up to 60% of the can- 
didates might not be confirmed. Another possibility is that 
most of the faint objects are not individual galaxies but un- 
resolved groups, a point made by Mori & Umemura (2006), 
who have shown that a number of small galaxies undergoing 
mergers can be identified as a single LAE at high redshifts. 
At 2 ~ 6.6, the resolution of the data is about 1", which cor- 
responds to a physical separation of 5.4 kpc. Although in our 
simulations, we do not find any objects that are separated 
by such small distances, this might also be due to insufficient 
resolution on scales of about 5 kpc. In this work, we discard 
this possibility, since, although a small number of such pairs 
might be found with higher resolution, they would probably 
not be enough to boost up the low luminosity end of the LF 
substantially. 

We now come back to clustering. The reason why its 
effects become important at z ~ 7.6 is easily explained us- 
ing Fig. [3] First, from the uppermost three panels (al, bl, 
cl), it is seen that the UVB photoionization rate (Tb) is 
very low at z ~ 7.6 (2.8 x 10~ 1B s" 1 ). Between z ~ 7.6 and 
6.6, it increases rapidly by about 3 orders of magnitude such 
that at 2 ~ 6.6, F B = 1.9 x 10" 13 s" 1 . Afterwards, the ion- 
ization rate increases only by a factor of about 2.5 between 
2 ~ 6.6 and 5.7. Second, at z ~ 5.7 6.6 the value of the 
boost parameter, (B) = 1+ (To) /F B [j is l ess than a factor 
of 1.5, while at z ~ 7.6, (B) w 58, i.e. the photoionization 
rate is strongly dominated by the local emission from the 
clustered LAEs. Since the IGM is already highly ionized at 
2 < 7 in the ERM, the extra local contribution from clus- 
tered LAEs does not affect the Lya transmission sensibly, 
as seen from the comparison of panels (a2,a3) and (b2,b3) 
of Fig. [3] However, at z ~ 7.6, the effects of clustering on 
the transmissivity, T a , are dramatic. The photoionization 
rate boost due to clustered LAEs makes the IGM transpar- 
ent enough that about 136 in our simulation volume become 
visible as compared to 30 that would be detected as LAEs in 
the absence of clustering effects. This is clear from the com- 
parison of panels (c2) and (c3), from which we conclude that 
the transmissivity is increased up to values of T a — 0.2 — 0.5 
when clustering is included in the computation as compared 
to T a = 0.2 — 0.3 when it is not. 

A few additional points are worth mentioning here: (a) 



3 (Tq) = (1/N) y\_ F(j(i), where rg(j) is the photoionization 
boost due to clustered LAEs seen by the i th LAE of the total N 
LAEs at the redshift considered. 
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Figure 3. Upper row: Photoionization rates as a function of LAE 
stellar mass, M,. The horizontal line corresponds to the contri- 
bution of the UVB, Tg; symbols in each panel denote the value 
of + Tq. Middle: Lyct transmissivity, T a including cluster- 
ing effects. Lower. Lya transmissivity without clustering effects. 
Columns refer to z ~ 5.7, 6.6, 7.6 as indicated. 

on a galaxy to galaxy basis, it is not necessarily the en- 
vironments of the most massive galaxies that experience 
the largest enhancement of the photoionization rate due 
to clustering as seen from the upper three panels of Fig. 
[3] This shows that the boost in photoionization rate seen 
by a small galaxy within the H n region of a more mas- 
sive one is greater than or comparable to the boost seen 
by a massive galaxy due to the contribution from a large 
number of small galaxies embedded in its H n region, (b) If 
Xhi <C 1, both including/excluding the contribution of Fg, 
T a increases with M 4 and hence with star formation rate 
(panels a2, a3, b2, b3). This leads us to conclude that for a 
highly ionized IGM, the contribution of T_l(oc M*) always 
dominates over that of Fq- (c) If Xhi w 1, T n still increases 
with the star formation rate both including/excluding Fg, 
although the scatter is much larger in the case including Fg- 
This is because excluding Fg, very few galaxies are luminous 
enough to transmit enough of the Lya luminosity to be visi- 
ble as LAEs (panel c3) . However, including Fg dramatically 
increases the T a such that about 4 times as many galaxies 
become visible (panel c2). Hence, we conclude that in this 
case, Fg contributes significantly in making the environment 
around LAEs more transparent to Lya photons. 

5.2 UV LF 

For each galaxy identified as a LAE included in the Lya 
LF, we calculate the total continuum luminosity in a band 
between 1250-1500 A, centered at 1375 A. We bin the num- 
ber of galaxies on the basis of the continuum luminosity 
magnitude and divide by the volume of the simulation box 
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Figure 4. UV LAE LF for the ERM. Points represent the data 
at two different redshifts: z ~ 5.7 (Shimasaku et al. 2006) (cir- 
cles), z nj 6.6 (Kashikawa et al. 2006) (triangles). Lines refer to 
model predictions including clustering at the redshifts (from top 
to bottom): z ~ 5.7,6.6, 7.6, for the parameter values in Tab. [l] 
The vertical dashed (dotted) lines represent the observational 2cr 
(3<r) limiting magnitude for z = 5.7 (z = 6.6). The shaded region 
in all panels shows the poissonian errors. Because radiative trans- 
fer is not accounted for in the simulations, at z = 7.6 our results 
represent lower limits to the observed objects. 

(75 3 /i~ 3 comoving Mpc 3 ) to obtain the UV LF for the LAEs 
identified in Sec. 15.11 This is shown in Fig. [4] 

We find that only a certain fraction of the continuum 
photons, f c , must escape the galaxy undamped by dust in 
the ISM to match to the observations; this fraction is found 
to decrease with decreasing redshift, going from 0.37 at z ~ 
6.6 to 0.22 at z ~ 5.7. This means that while 37% of all 
continuum photons escape the galaxy at z ~ 6.6, only 22% 
escape at z ~ 5.7. We interpret this decreasing f c to be the 
result of an increase in the dust content of the galaxies. To 
make predictions at z ~ 7.6, we again use the same value 
of f c as at z ~ 6.6, making the assumption that the dust 
content of LAEs does not evolve between these two redshifts. 
The values of f c for all the three redshifts are shown in Tab. 

m 

An interesting point here is that at z ~ (5.7,6.6), 
fa/ fc ~ (1.4,0.8), confirming a result obtained from our 
previous semi-analytical model, presented in DFG08. Using 
the Supernova extinction curve, f a /fc ~ 0.8 for a homo- 
geneous distribution of dust, Bianchi & Schneider (2007). 
However, no single extinction curve (Galactic, Small Mag- 
ellanic Cloud or Supernova) can give a value of fa/fc > 1- 
The relative damping at z ~ 5.7 can only be explained by 
an inhomogeneous two-phase ISM model in which clumped 
dust is embedded in a hot ISM, as proposed by Neufeld 
(1991). The data thus, seem to hint at the fact that the 
dust distribution in the ISM of LAEs becomes progressively 
inhomogeneous/clumped with decreasing redshift. 

We calculate the color excess E(B — V) (= A v /R v ) for 
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Figure 5. Comparison of theoretical SEDs (lines) with observations (points) for LAE #07, #08 and #34, (from left to right) from Lai 
et al. (2007). Points with downward pointing arrows represent the 3<r upper limits of the data. 



each of the redshifts as a sanity check. Using the supernova 
dust extinction curve (Bianchi & Schneider 2007), we calcu- 
late A v = A A (1375 A)/5.38 and R v = 2.06. The color excess 
is then calculated as 



E(B-V) 



A A (1375A) 
5.38R V 



2-51og 10 f c 
5.38R V 



(9) 



With the above formulation, we find that E(B — V) ~ 0.15 
atz~ 5.7, E(B - V) ~ 0.1 at z ~ 6.6. These finding are 
consistent with other independent data as well: in fact, Lai 
et al. (2007) have shown that fitting the SEDs of 3 LAEs ob- 
served by them at z ~ 5.7 requires E(B-V) < 0.425 (0.225) 
for young (old) galaxies. 



Table 2. Physical properties of the LAEs from the simulation 
which match the observed Lya luminosities most closely. For each 
LAE observed by Lai et al. (2007) (col. 1), we show the age (col. 
2), stellar metallicity (col. 3), SFR (col. 4) and color excess (col. 
5) for the corresponding LAE from our simulation. 



#LAE 


U 


Z 


M, 


E(B - V) 




(Myr) 


(20 ) 


(M yr-l) 




07 


191 


0.23 


9.7 


0.15 


08 


182 


0.32 


9.6 


0.15 


34 


220 


0.23 


7.3 


0.15 



5.3 SEDs 

Lai et al. (2007) have observed the spectra for three LAEs 
(#07, #08 and #34) at z ~ 5.7 with observed Lya lumi- 
nosities, L a = (4.9,4.3,3.6) x 10 42 ergs _1 (K. Lai, private 
communication). From amongst the galaxies we identify as 
LAEs at 2 ~ 5.7, we select those three galaxies whose ob- 
served Lya luminosities match most closely with the three 
values observed. The spectra of each of the selected galax- 
ies (obtained using STARBURST99, as explained in Sec. [3j) 
is then attenuated using the SN dust extinction curve and 
E{B — V) = 0.15, to produce synthetic spectra which are 
shown in Fig. [5] 

As shown by Lai et al. (2007), the spectra can be well 
fit by different kinds of stellar populations, varying the age, 
metallicity and color excess. However, we have no free pa- 
rameters since the age, metallicity and SFR for each of the 
galaxies is obtained from the simulation outputs and the 
color excess value is obtained using the UV LF at z ~ 5.7 
(Shimasaku et al. 2006) as mentioned before. As seen, a 
remarkable agreement is found between the synthetic and 
observed SEDs. Even though we are using a small and bi- 
ased data set of LAEs selected based on the IRAC 3.6/ira 
and 4.5/xm detections, the agreement between the synthetic 
and observed SEDs provides a strong consistency test of our 
model. 

We briefly mention the physical properties of the LAEs 
whose synthetic spectra are shown in Fig. [5] The stellar 
ages are 182-220 Myr for the three LAEs (see also Tab. 



[2|. Hence, these LAEs are intermediate age objects, rather 
than being very old (t* ~ 700 Myr) or very young (t» ~ 5 
Myr). Their stellar metallicities are about 0.2 — 0.3 Z@, and 
the SFR are between 7-10 M Q yr _1 . Further, the SEDs for 
all the three LAEs are well reproduced by a single value of 
the color excess which shows that all these objects possibly 
contain similar amounts of dust in the ISM. 



5.4 The nature of LAEs 

Having selected the galaxies in the simulated volume that 
would be experimentally defined as LAEs, and having shown 
that their luminosity function and SED are consistent with 
data, we can take a further step and quantify the physi- 
cal properties of these objects. This includes the relation 
between stellar and halo mass, age, metallicity, star forma- 
tion rate and evolution. We will present these results up to 
z ~ 7.6; however it must be kept in mind that at that red- 
shift, only 136 galaxies are bright enough to be observed us- 
ing standard narrow-band techniques which makes the scat- 
ter much larger as compared to the lower redshifts. 

The first set of relations among stellar and halo mass, 
age, metallicity is presented in Fig. 6. LAEs are character- 
ized by dark matter halo masses in the range 1O 1O ' 2_11,8 M0 
at = 5.7, corresponding to > 2<r fluctuations at all red- 
shifts. This range becomes progressively narrower at earlier 
epochs because of two occurrences: (i) the mass function 
of halos shifts to lower masses in hierarchical structure for- 
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Table 3. For all the LAEs comprising the Lyct LF including clustering effects, at the redshifts shown (col. 1), we show the range of mass 
weighted ages (col. 2), the average stellar metallicity (col. 3), the average SFR (col. 4), the average intrinsic and observed EWs (col. 5, 
6), the average value of the star formation indicator (col. 7) and the SFR density (col. 8) obtained from our simulation. 
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Figure 6. Physical properties of LAEs at z ~ 5.7, z ~ 6.6 and 
z ~ 7.6 as referred to by columns. As a function of stellar mass, 
M„, we show (a) the halo mass, M h , (al-cl), (b) mass- weighted 
stellar ages, i„, (a2-c2), and (c) mass- weighted stellar metallicity, 
Z, (a3-c3). 



mation models; and (ii) smaller halos progressively become 
invisible at higher redshifts as their stellar mass and hence 
luminosity is too low to be detected (we recall that all the 
objects analyzed are part of the LFs shown in Figs. [2] and [4]). 
In spite of this fact, there is almost no evidence of evolution 
of the stellar to halo relation in the redshift range under ex- 
amination, and we find M* oc M^' 56 , with the best fit given 
by the following expression: 

log 10 (A/ h ) = (0.64-0.06Az) log 10 (M*) + (5.0+0.5Az), (10) 

where Az = (z — 5.7). The above relation implies that the 
stellar mass per unit halo mass is increasing towards larger 
systems, which is just a restatement of the well-known fact 
that the star formation is less efficient in small galaxies due 
to the inhibiting effects of mechanical feedback. The typical 
stellar masses of LAEs are < 10 10 ' 5 at all redshifts, i.e. they 
are smaller than the Milky Way. 

Although a considerable spread is present, on average 
larger systems tend to be older than small ones. LAEs re- 
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Figure 7. (a) SFR, M», and (b) SFR indicator, X = M,t»/M», 
as a function of stellar mass (M*) for LAEs at three redshifts: 
z ~ 5.7 (circles), z ~ 6.6 (triangles) and z ~ 7.6 (asterisks). 



suit from the merging of several sub-halo systems in which 
star formation was ignited long before. In Tab. [3] we show 
that the range of ages reduces towards high redshifts from 
38 - 326 Myr at z ~ 5.7 to 39 - 165 Myr atz~ 7.6. These 
ages imply that the oldest stars in these systems formed al- 
ready at 2 = 8.5, thus during the reionization epoch. We 
reiterate that, because of finite numerical resolution, stars 
in halos below our resolution limit might have formed even 
before that epoch, so the previous value must be seen as 
a lower limit to the onset of star formation activity. From 
the previous discussion it is clear that LAEs are neither 
pristine nor very young objects forming their stars for the 
first time. This is confirmed by the range of metallicities 
found in these systems and shown in Fig. [6] (see also Tab. [3]). 
Metallicities are in the range 0.02 — 0.55 Z@, with the mean 
over the sample decreasing with redshift; also, larger objects 
have higher metallicities as expected from their larger stellar 
masses. These results are broadly consistent with an enrich- 
ment predominantly caused by SNII, and therefore tightly 
following the SFR of the galaxy; the increasing scatter seen 
at the lowest redshift hints at a larger contribution by SNIa. 
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Figure 8. Contour plot showing relation between observed Lyc* luminosity, age and SFR at z ~ 5.7. The SFR values, Log(M* /Mq yr x ), 
are color-coded by the bar on the right. 



The mass-metallicity relation for LAEs is also an interesting 
outcome of our study. This can be conveniently expressed as 

Z/Z Q = (0.25 - O.O5A2) log 10 (M») - (2.0 - 0.3Az), (11) 

where Az = (z — 5.7). Differently to the analogous relation 
observed at lower redshifts (Tremonti et al. 2004; Panter 
et al. 2008; Maiolino et al. 2008), we do not see the sign 
of a flattening of metallicity towards larger masses. As the 
flattening is usually interpreted as a result of a transition 
from a galactic wind-regulated metal budget to a close-box 
evolution in which all metals are retained, we conclude that 
in our LAEs, winds play an important role because of their 
relatively low masses, a point already noticed above. 

We now turn to the analysis of the star formation prop- 
erties of LAEs (Fig. [7}. Star formation rates are in the range 
from 2.5-120 Moyr -1 , i.e. a sustained but not exception- 
ally large star formation activity, with large objects at later 
times being the most prominent star factories. The relation 
appears to flatten below M* ^ 10 9,3 Mq, as is clearly seen. 
Stated differently, low mass LAEs have star formation rates 
confined in the narrow range 2.5-10 Mq yr -1 , whereas only 
(relatively few) larger objects undergo intense star formation 
events, M, > 50 Mq yr -1 . This might be the result of feed- 
back regulation, which prevents small objects from burning 
most of their gas fuel at high rates. To investigate this aspect 
more, we have studied the behavior of a star formation in- 
dicator, I = M*t*/M*. Physically, this is the ratio between 
the stellar mass produced if a LAE had always formed stars 
at the rate deduced at the given redshift for a time equal to 
the mean age of its stars, and the actual total stellar mass of 



the system. Hence if 1 > 1 (I < 1), the star formation rate 
was lower (higher) in the past. The large majority of LAEs 
show values of 1 < 1, indicating that the star formation rate 
averaged over the entire history must be larger than the final 
value. This could be due to huge bursts of star formation, 
catalyzed by mergers. Some systems (at the low and high 
ends of the stellar mass distribution at z ~ 5.7) have I > 1, 
i.e, they have quietly built up their stellar population at an 
increasing rate as they grow in mass by subsequent merg- 
ings. We provide a handy fit for the SFR in terms of the 
stellar mass: 



Iog 10 (M.) = (0.4-0.04A 2 )[log 10 (M*)] 2 - 
(6.9 - O.SAz) log 10 (M„) + (30 



5Az) 



where Az = (z — 5.7). 

Finally, we discuss the dependence of the observed Lya 
luminosity on age and SFR (Fig. [8}. The SFR for all galaxies 
with a given L a are very similar, even though the ages vary 
between 40 — 340 Myr. As explained before, this is because 
according to the ERM, at z ~ 5.7, the Universe is already 
so ionized (xbi ~ 10~ 5 ) that even the smallest emitters 
are able to build large enough Stromgren spheres on short 
timescales. Further, the number of galaxies with small SFR 
(M* ^ 20Moyr -1 ) is quite high, after which this number 
decreases rapidly (see also Fig. [7}, with very few galaxies 
in the highest luminosity bins. Thus, the faint end of the 
LF samples a mixture of young and old objects, while LAEs 
in the bright end are predominantly massive, intermediate 
(200-250 Myr) age systems. 
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Figure 9. Normalized distribution of number of LAEs as a func- 
tion of the mass weighted stellar age (t*). The lines are for z: 5.7 
(solid), 6.6 (dotted) and 7.6 (dashed). 

6 DISCUSSION 

The ages of LAE stellar populations are currently hotly de- 
bated, hence we need to discuss our results in the frame- 
work of the various arguments given in the literature. For 
example, Finkelstein et al. (2009) find a bimodality in their 
sample of 14 LAEs at z ~ 4.5; their objects are either very 
young (< 15 Myr) or very old (> 400 Myr). Their results 
could be explained by invoking two different star formation 
modes in LAEs: a recent strong burst or a continuous SF 
in which the bulk of the population is dominated by old 
objects. To better evaluate this possibility, we plot the nor- 
malized distribution of the mass- weighted stellar ages for all 
the LAEs identified in the simulation volume at z ~ 5.7, 6.6 
and 7.6 in Fig. [§] We do not find any such bimodality from 
our simulation. LAEs are distributed in age between 39-165 
Myr at z ~ 7.6 and this range increases with decreasing red- 
shift, as already mentioned above. A possible explanation 
of the experimental result could be the presence of dust. 
LAEs would be visible in the Lya when the age < 15 Myr 
so that not enough dust would have formed and at later 
times > 400 Myr when the galaxy would have destroyed/ 
removed most of its dust content. In the intermediate pe- 
riods, the Lya line would be attenuated below observable 
limits. However, a complete picture is possible only if dust 
creation/destruction is modelled accurately, which we leave 
for future work. 

Another interesting point is the Lya line equivalent 
width. We calculate the rest frame observed EWs at z ~ 5.7 
using our estimates of f c = 0.22 and f a — 0.3 from the UV 
and Lya LFs as mentioned in Tab. [1] as 

EW = EW mt ^(l- f esc )T a , 

Jc 

where EW int = L™* / L™* (1375 A). For our sample we find 
that 55 < (EW/A) < 89, with mean value of 63. 5A. This 
is about a factor of 2 lower than the mean value of 120 A, 
observed by Shimasaku et al. (2006). This discrepancy is 
probably due to physical effects that are not included in our 
model and which could broaden the EW distribution and 
increase its mean value. These include: (i) Lya emission 
from cooling of collisionally excited H i in the ISM (Dijk- 
stra 2009). Adding this contribution to the Lya luminosity 
from stars could increase the intrinsic Lya luminosity from 
the galaxy by large amounts while leaving the intrinsic con- 



tinuum luminosity value unchanged, thereby increasing the 
intrinsic EW; (ii) we are using the same escape fractions of 
Lya and continuum photons for all the galaxies. In reality, 
however, these values would change on a galaxy to galaxy 
basis, depending on the amount, topology and distribution 
of dust inside the galaxy, contributing to a broadening of the 
EW distribution; (iii) outflows could increase the observed 
EW distribution since they would enable more of the Lya 
luminosity to escape undamped by dust, while leaving the 
continuum unchanged. Further, outflows with a sufficiently 
high value of the H i column density (~ 10 18 " 21 cm" 2 ) add 
a bump to the red part of the Lya line, thereby increasing 
the observed EW ( Verhamme et al. 2006) ; (iv) inflows would 
lead to a further damping of the Lya luminosity and lead 
to a decrease in the EW. A full exploration of these effects 
requires modelling the amounts and distribution of dust and 
peculiar gas motions in individual galaxies. We defer this to 
future works at this point. 

At 2 ~ 5.7, we find that for all galaxies with an ob- 
served Lya luminosity, L a ^ 10 42,2 ergs _1 , the observed 
EW is greater than 20A. This means that all LBGs having 
a Lya luminosity larger than this value show an observed 
Lya EW > 20 A from our work. Although their samples 
are very small (6 and 8 galaxies respectively), Stanway et 
al. (2004) and Dow-Hygehmd et al. (2007) find that only 
about 30% of LBGs showing a Lya emission line have an 
observed Lya EW larger then 20 A. The main reason for 
this tension could be the fact that we are using the same 
ratio of fa/fc = 1-4 for all the galaxies at z ~ 5.7; how- 
ever, this value should change on a galaxy to galaxy basis. 
The relative escape fractions of Lya and continuum pho- 
tons must depend on the amount and distribution of dust 
and H i inside the galaxy. Modelling this effect more realisti- 
cally would probably lead to many galaxies having EW < 20 
A and resolve this tension. 

LAEs are among the most distant galaxies known. As 
the quest for reionization sources is struggling to iden- 
tify the most important populations for this process, it 
is worthwhile to assess to what extent LAEs might rep- 
resent such a long-searched-for population. First we note 
that (Tab. 0) at z ~ 5.7, the SFR density provided by 
LAEs is p, = 1.1 x 10" 2 M yr" 1 Mpc" 3 , and this de- 
creases with increasing redshift. Comparing this to the cos- 
mic SFR density of 0.12 M yr" 1 Mpc -3 at the same red- 
shift measured by Hopkins (2004), we find that LAEs con- 
tribute only about 9.2% to the cosmic SFR value, confirm- 
ing the previous result of DFG08. Corresponding to this 
SFR evolution and / esc = 0.02, it is easy to derive that 
the H i ionizing photon rate density (qlae) contributed by 
LAEs is q LA E = 3.1 x 10 49 s" 1 Mpc" 3 at z = 6.6. We com- 
pare this photon rate density to the one necessary to balance 
recombinations given by Madau, Haardt and Rees (1999): 

where C is the IGM clumping factor. For (7=1, correspond- 
ing to a homogeneous IGM, we get the minimum H i ionizing 
photon rate density necessary to balance recombinations, i.e. 
q rec = 1.24 x 10 50 s" 1 Mpc" 3 . This means that the LAEs on 
the Lya LF at z ~ 6.6 can contribute at most 25% of the 
H i ionizing photons needed to balance recombinations at 
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this redshift. This value decreases to 0.8% as C increases to 
30. 

The total H i ionizing photon rate density contributed 
by all astrophysical sources (q a u) at z ~ 6.6 can be calcu- 
lated as (Bolton & Haehnelt, 2007) 

g o „ = 10 5113 r- 12 (^^^ s^Mpc" 3 ; (13) 

using T_i2 = 0.19 (Tab. [[]), gives q all = 2.5 xlO 50 s" 1 Mpc -3 
at z ~ 6.6. This implies LAEs contribute about 12.5% of 
the total H i ionizing photon rate, which is consistent with 
the SFR density contribution of LAEs to the global value. 
Such a low contribution from LAE, given their halo and 
stellar masses does not come as a surprise. Indeed, this result 
is consistent with the previous estimates by Choudhury & 
Ferrara (2007) who showed that only a fraction £ 1% of 
the photons required to ionize the IGM come from objects 
in the LAE range, the bulk being provided at high redshifts 
by faint (or even Ultra Faint, see Salvadori & Ferrara 2009) 
dwarf galaxies. 



7 SUMMARY 

We use a large scale hydrodynamical simulation to derive the 
physical properties (halo, stellar, gas mass, star formation 
rate, stellar age and metallicity) of high redshift galaxies. 
Using a reionization history consistent with an early reion- 
ization epoch (ERM), we fix the evolution of the background 
photoionization rate, and calculate the size of the Stromgren 
sphere around each object, taking into account the ioniza- 
tion boost due to spatial clustering of the sources when com- 
puting the corresponding Lya radiation transmission. We 
then define as LAEs galaxies with an observed Lya luminos- 
ity above the currently observable limit (L a ^ 10 42,2 erg s _1 ) 
and an observed EW > 20 A. 

We find that Lya photons at both z ~ 5.7 and 6.6 must 
be attenuated by dust within the galaxy to reproduce the 
observed LAE luminosity function. The amount of photons 
that escape the ISM, undamped by dust is f a = 0.3 at 
both redshifts. At z ~ 6.6 the simulated LF lies between 
the upper and lower limits of the data. At face value, this 
results implies that further spectroscopy of LAE candidates 
at this redshift should rule out ~ 60% of them as LAEs. 

Clustering of sources boosts the average value of the 
total photoionization rate (and consequently the ionized hy- 
drogen fraction) in the surroundings of galaxies by < 1.5 
times at z ~ 5.7, 6.6; this value increases by more than 
a factor of 50 at z ~ 7.6. For an almost neutral IGM 
(z ~ 7.6, xhi = 0.16), if clustering effects are ignored, only 
those few galaxies that are able to carve out a large enough 
H ii region are visible. However, about 4 times more galax- 
ies become visible when clustering effects are included. In a 
highly ionized IGM (z ~ 6.6,5.7, xhi ~ 10- 4 ,10~ 5 ), how- 
ever, the effect of clustering is indiscernible on the Lya LF. 
Note that, the clustering boost is not necessarily the highest 
for the most massive galaxies. Our model is probably some- 
what overstating the effects of clustering because we are not 
treating the radiation transfer of ionizing radiation. 

The fraction of Lya luminosity transmitted always in- 
creases with SFR and about 50% of the Lya luminosity is 
transmitted at z = 5.6, 6.6; however, this relation shows a 



larger scatter at z — 7.6 (transmission between 20-50 %) 
compared to the lower redshifts since at z ~ 7.6, it is es- 
sentially the photoionization rate boost due to clustering 
which enables Lya luminosity transmission, if galaxies are 
embedded in the neutral patches of the IGM. 

The presence of dust even at these high redshift is fur- 
ther confirmed by the analysis of the LAE UV luminosity 
function, which we have also reproduced within our study. 
Matching the UV LF requires that the escape fraction of 
UV continuum photons is f c — (0.22, 0.37) at z = (5.7, 6.6), 
again a clear signature of dust attenuation. We interpret 
the higher attenuation of the continuum photons relative 
to the Lya (Neufeld 1991) at z ~ 5.7 as hinting at an in- 
homogeneous two-phase ISM with dust clumps embedded 
in a warm intercloud gas while the relative attenuation is 
consistent with homogeneously distributed dust at z ~ 6.6. 
At z ~ 5.7, the color excess calculated from our model, 
E(B — V) — 0.15, is consistent with the value deduced by 
Lai et al. (2007), E(B-V) < 0.425, and the results obtained 
by Nagamine et al. (2007) (E(B - V) = 0.15) using SPH 
simulations. The color excess, E(B — V) = 0.1, at z ~ 6.6 is 
lower than the value at z ~ 5.7; this is expected due to to the 
average lower metallicities of LAEs towards higher redshifts 
(see Tab. 3). By extrapolating these results at z = 7.6, we 
have predicted the UV and Lya LFs for LAEs at z ~ 7.6. 

As a further test of our results, we have selected three 
LAEs from our simulation whose observed Lya luminosities 
match most closely with the values observed at z ~ 5.7 by 
Lai et al. (2007). For each of the chosen objects, the intrin- 
sic spectra (obtained using STARBURST99) is attenuated 
using a SN dust extinction curve and E(B — V) = 0.15 to 
match to the observations. We find our synthetic Spectral 
Energy Distributions are in very good agreement with the 
observations. The objects observed are found to have ages 
of a few hundred Myr, metallicity about 20 — 30% Zq and 
SFR of 7-10 M yr" 1 . All three LAEs are well fit by a single 
value of the color excess, hinting at the fact they all might 
be dust enriched to a similar degree. 

Finally, we have discussed in detail the physical prop- 
erties (and we provide handy fitting functions to several re- 
lations among them) of galaxies which we identify as LAEs 
at z = 5.6, 6.6 and 7.6. The ages of the LAEs range between 
39 — 165 Myr at z ~ 7.6 but this range increases to 38 — 326 
Myr at z ~ 5.7. Further, the average metallicity of LAEs 
is Z = 0.12^0 at z ~ 7.6, and increases to Z = 0.22Zq at 
z ~ 5.7. Hence, LAEs are more metal-enriched than what 
usually assumed. Star formation is relatively suppressed by 
feedback in low mass halos (M h < 10 n M o ) and it rises 
steeply for the larger halo masses. The large majority of 
LAEs seem to have had a higher average star formation rate 
over the entire history as compared to the final value. Fur- 
ther study would be necessary to put this statement on more 
solid grounds. 
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